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Supersymmetric loop corrections induce potentially large CP- violating couplings of the Higgs bosons 
to nucleons and electrons that do not vanish in the limit of heavy superpartners. The Higgs-mediated 
CP-odd four-fermion operators are enhanced by tan 3 /3 and induce electric dipole moments of heavy 
atoms which exceed the current experimental bounds for the electroweak scale Higgs masses and 
tan /3 > 10. If only the first two sfermion generations are heavy, the Higgs-mediated contributions 
typically dominate over the Barr-Zee type two-loop diagrams at tan f3 > 30. 
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Non-observation of the electric dipole moments 
(EDMs) of the neutrons yj] and heavy atoms poses 
a serious problem for low-energy supersymmetry. This 
is because generically the EDMs are induced already at 
the one loop level and their predicted values exceed the 
experimental limits by orders of magnitude (see || for 
recent analyses). The most straightforward way to cir- 
cumvent this problem is to assume that the superpartner 
mass scale is rather high (over a few TeV) which leads to 
the suppression of all effective one-loop-generated CP- 
odd operators of dimension five and higher, and seem- 
ingly allows for arbitrary CP-violating phases. The aris- 
ing Higgs mass fine-tuning problem can be alleviated if 
the third generation sfermions are kept light. 

In this work, we consider sources of the EDMs which 
survive the decoupling of the superpartners while the 
Higgs masses are kept fixed. In this limit, the Minimal 
Supersymmetric Standard Model (MSSM) degenerates 
into a two-Higgs doublet model (2HDM) with an impor- 
tant one-loop-induced modification of the Yukawa sector 
compared to the usual type II models Q . We show that 
the Higgs-mediated four-fermion operators induce poten- 
tially large electric dipole moments of heavy atoms that 
grow as tan 3 (3, and that supersymmetric models with 
large tan /3 face the SUSY CP problem even in the limit 
of heavy SUSY particles. We further compare the Higgs- 
mediated contributions with the two-loop induced con- 
tributions |?],H which are significant if the masses of the 
third generation squarks are near the electroweak scale. 

At the tree level, the down type quarks and charged 
leptons obtain their masses from the interaction with the 
first Higgs doublet Hi . The finite one loop SUSY correc- 
tions induce considerable couplings of the second Higgs 
doublet H2 to the D-quarks and charged leptons, that 
are absent in the limit of unbroken supersymmetry: 

-C Y = Y^HxDlDr + y D h\d l d r 
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where Y^ E are the tree level Yukawa couplings and ~S>d,e 
are the loop induced couplings. The typical representa- 
tives of important threshold diagrams are given in Figs, la 
and lb. Ignoring flavor changing effects, 
and y D E by: 



The loop functions Jr> and Je are complex and contain 
the dependence on the CP phases in the soft-breaking 
sector, which leads to the CP-odd interactions of the 
physical Higgses with the -D-quarks and leptons. When 
one redefines the phase of the right-handed D and E 
fields, D R -> e~ i5 °D R and E R -> e~ l5E E R such that the 
mass term in Eq. (Q) becomes real, the induced CP-phase 
in the Higgs interaction gets enhanced by tan/3 = V2/V1: 
Sd,e = Arg(l + Ju^tan/3). Thus, large tan/3 can com- 
pensate the loop smallness of Jd,e so that the phases 
5d,e can be order one. An exchange by physical Hig- 
gses will then produce CP-odd four-fermion interactions, 
Fig. 2a. The relevant interactions are induced by the ex- 
change of a CP-odd Higgs boson A and CP-even Higgs 
boson H between the D-quarks and the electron, and 
between the D-quarks: 
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SM Yf M (suxSi- sin Sj) 
Ji tan /3| 1 1 + Jj tan (3\ 

xipi-ipi i>jij 5 4>j . (3) 



Here Yf denote the Standard Model values for the 

Yukawa couplings, Yj? M = y^mj/v (v=246 GeV), and 
the CP-phases are understood modulo tt. In the deriva- 
tion of (||), we have used the relations ttia — ran, 
cos 2 a ~ 1, where a is the neutral Higgs mixing angle, 
and dropped all tan /3-suppressed terms. This approxi- 
mation works well even for moderately large tan (3. 



<H,> 



,<H,> 



A ; 



D 



(a) 



b : 



(b) 



\t,b 



\ H 



H 



(C) 



FIG. 1. SUSY threshold corrections in the down quark 
Yukawa sector, (a) and (b), and in the Higgs sector, (c). 

CP-odd contact interactions can also be induced via 
A — H mixing, Fig. 2b, which appears due to CP- 
violating Higgs couplings to the third generation squarks 
§ (Fig.lc): 
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y E = j e y e 0) 
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Here we have used m 2 H ~ m 2 A » (AH). Such effects 
were studied previously in the context of 2HDMs with 
spontaneous breaking of CP by Barr [jlO| . 

Inspection of Eq.(||) reveals that the CP-odd coupling- 
grows as tan 3 /3 because sin Si ~Im Ji tan /3/| 1 + Ji tan /3| , 
until the radiative corrections become comparable to the 
tree-level values. Here we treat to a as an independent 
variable (it is proportional to the SUSY Bfi parameter). 
The cubic growth is different from a tan 2 /3-behavior in 
2HDMs with spontaneous CP violation [fuJ]. Thus, gen- 
erally Eq.(^) represents a subleading effect, as (AH) con- 
tains a loop smallness not compensated by large tan/3. 
We note that the QCD renormalization group flow for 
the electron-quark interactions from niA to lGeV is triv- 
ial at one loop and one can simply take Yf M normalized 
at lGeV. 
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FIG. 2. Higgs-mediated four-fermion interactions with CP 
violation in the Higgs-fermion vertex (a) and on the Higgs 
line (b). 

Using Eq.(||), we calculate the EDMs of paramagnetic 
atoms and estimate the EDMs of diamagnetic atoms and 
neutrons. The semi-leptonic operators in (Q) induce two 
types of T-odd nucleon-electron interaction 



-CP 



C s NNeij 5 e + C P Ni-f 5 Nee 



(5) 



with possible isospin dependence. Cs and Cp are 
severely constrained by the recently improved experimen- 
tal bounds on the EDM of the thallium and mercury 
atoms |,|. With the use of the standard technique for 
the QCD matrix elements [jllj of a heavy quark over a 
nucleon state, the isospin-singlet coupling Cs can be ex- 
pressed as 

c ^ 5.5 x 10- 10 tan 2 ^ 
to 2 1 1 -\- J e tan /3| 



(1 — 0.25k) (sin 5b — sin<5 e ) 



3. 3k (sin S s — sin<5 e ) 
|1 + J s tan^| 



|1 + J 6 tan/3| 
0.5 (sin Sd — sin 
|1 + J d ta,n/3\ 



S e ) 
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Here we have used (m u + md)(N\uu + dd\N)/2 = 45 
MeV and (m„ - m d )(N\uu - dd\N)/§Q MeV< 1. The 
coefficient k = (N\m s ss\N) /220 MeV parametrizes the 
uncertainty in the value of (N\m s ~ss\N) matrix element. 
Its "best" value k = 1 is inferred from the leading order 
flavour SU(3) analysis of the baryon octet mass split- 
tings. An assumption that the strange quark behaves as 



a heavy quark would lead to a smaller value, k = 0.3. It 
is important to note that a significant source of uncer- 
tainty - the poorly known masses of the light quarks - 
does not affect Eq.([|)! Using (||), and the results of the 
atomic calculation that relates rf/ri and Cs [fl2|, 



2 T 1 



1.5 x 10 _17 e cm x C s (100 GeV) 2 



(7) 



one can express the thallium EDM in terms of the SUSY 
parameters. Comparison of (f^) with the experimental 
data provides the bound C s < 1.1 x 10~ 8 (100 GeV) -2 . 

The dimensionless loop functions JVs depend on the 
pattern of the soft masses. To get an idea of the size 
of the induced EDMs, let us first consider a toy model 
with TOsfcrmion = TO gaug ino = = | Ai \ = M > Mz- The 
dominant contribution comes from the squark-gluino and 
stop-Higgsino exchange: 



3tt 



J e — 0; Jd — J t 

(X 

Jb= tt~ exp{i0 M + i<j> 3 } 
on 



expji^ + i(j> 3 }; 



(8) 



1 2 ex P{^M + #A t }, 
OZ7T 



where 



, 4>A t are the phases of the /^-parameter, the 



gluino mass, and the At parameter, respectively. In the 
case of general soft terms, the gluino contribution to J L 
should be multiplied by |/iM 3 |/(TO 2 1 , to 2 2 , |M 3 | 2 ) and the 
Higgsino contribution by \^A t \ /(m 2 i, to 2 2 , |/i| 2 ) 3], where 
mu 2 are the squark mass eigenvalues, M 3 is the gluino 
mass, and the loop function is defined by 

ablnialb) + bcAnib/c) + ac\n(cla) 

I(a,b,c) = 2 — (9 

(a — b)(b — c)(a — c) 



such that / = 1/A/ 2 for a = b 



M 2 . In the same 



limit, the CP-odd Higgs mixing is given by 



(AH) = 



SMn4 c 



sin(2<£ M + 2<t> At ) 



64tt 2 

+ (V b 5M ) 4 tan 4 /3sin(2^ + 2$ A J] 



(10) 



Obviously, both Jj's and (AH) are independent of the 
superpartner mass scale M. An expression for (AH) in a 
more general case can be found in ^ . As we will see, the 
effect of the A — H mixing does not impose significant 
constraints for rriA > 150 GeV, so henceforth we will 
mainly concentrate on the effect of the vertex corrections. 

The Higgs-quark vertex corrections lead to the follow- 
ing thallium EDM normalized to the current 90% C.L. 
experimental bound [e?Ti]cxp = 9.4 x 10~ 25 e cm j2j: 



d T i 



tan 3 /? 



Mti]cx P 350 x m 2 00 



[sinc/> M +0.04sm(</> M + </> At )], (11) 



where we have set 3 = 0, |1 + J;tan/3| ~ 1, and k = 1. 
TOioo is 771,4 measured in the units of 100 GeV. Already at 
tan/3 ~ 7 the r.h.s. of ( |Ti"| ) may reach 1 (while the stan- 
dard EDM contributions are suppressed for the SUSY 
masses of 10 TeV). 
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We conclude that even for arbitrarily heavy su- 
perpartners, the SUSY CP problem reappears if 
tan/3(mA/100 GeV) 2 ^ 3 > 10. For instance, with </y ~ 
1, m A ~100 GeV, and tan ,3 ~ 60 the induced EDM ex- 
ceeds the experimental bound by almost three orders of 
magnitude! It is important to note that these calcula- 
tions are free of large nuclear uncertainties J^] . 

The CP-odd constant Cs also induces EDMs of dia- 
magnetic atoms through the mixing with the hyperfine 
interaction 13 Our prediction for the mercury EDM 
is 



^Hg(Cs) 
[^Hgjcxp 



tan 3 /? 
900 x m 2 00 



[sin^ + 0.04sin(^ + ^A t )], (12) 



where the current experimental bound is [c?Hg]oxp = 
2 x 10~ 28 e cm ||. This imposes a slightly weaker bound 
that the thallium EDM. However, in the case of diamag- 
netic atoms there are two additional classes of contribu- 
tions, induced by Cp and the nuclear Schiff moment. To 
evaluate Cp, we follow the strategy of Ref. |p^| . Unlike 
the previous case, there is a strong dependence of the re- 
sult on m u /rrid, and within the error bars for this ratio 
the matrix elements of 5175 s and 6175 b over the neutron 
are compatible with zero. The d-quark contribution gives 



Cp(n) 



6.3 x 10" 



tan 2 /3(sin 84 — sin S e 



|1 + J e tan/3||1 + J d tan/3|' 



(13) 



where we have used m u /rru = 0.55. Using the results 
of the atomic calculation [Q, we convert this into the 
bound 



[cfeglcx 



tan 3 /? 



3500 x m 2 00 



(14) 



This is clearly a subleading contribution, compared to 
(|l2l). Finally, the 5*759 99 interactions in (||) induce c?Hg 
via the nuclear Schiff moment. It is known that a T- 
odd one pion exchange between nucleons is the domi- 
nant source of the Schiff moment. We estimate the T- 
odd pion-nucleon coupling g^NN following the approach 
of |g5[. Since sinS s — smSj in our case, there are can- 
cellations in the sum (^) and the dominant contribution 
comes from the operator dij^d bb. The resulting isospin- 
triplet coupling g^ NN ir NN is 



9ttNN 



2.7 x 10 



_ 13 tan 2 (su\8d — sind{,)(l — 0.25k) 
m ioo |1 + Jbtan/3||1 + Jdtan/3| 



Skipping a long chain of nuclear and atomic matrix ele- 
ments that relate g^^N an d ^Hg (see Ref. [|l2| for details), 
we get: 



rf Hg (Schiff) 

cxp 



tan 3 /3 



1.1 x 10 4 m 2 00 



sin(^ + i/>A t ). (15) 



Remarkably, Eq.(|15[) has the same sensitivity to </> M + 
(f>A t as Eq. ( [ill ) > y et this calculation involves considerable 



uncertainties. A combination of dxi and dn g constrains 
both phases, (f>^ and <f>A t , once again exemplifying the 
complementarity of the two measurements Jl6| . 

Let us now consider a popular scenario where only the 
first two sfermion generations are assumed to be heavy 
(> 10 TeV). We further assume the most conservative 
case when the gluinos are also heavy such that the leading 
term in Eq.(|ll|) disappears. Then we have 

, 2 



Jd — J s — J e — 0; Jb 



{Y SMy 



A t n I(m 2 ii7 mf 2 ,\fi\ 2 ) 



32tt 2 

We note, however, that in order to suppress the contri- 
bution of the gluino exchange diagram in Jb one would 
have to require M3 > 60^! 

Fig. 3 shows a tan/3 dependence of Cs/[Cs] cxp = 
dxi/[^Tl]exp in this model. We use the parameters of Ref. 
g, namely \A t , b \ = \fi\ = 1 TeV, Arg(^ t , b ^) = tt/2, and 
77i sq = 0.6 TeV (assuming that the left and right stop 
and sbottom mass parameters are given by m sq while 
other squarks are decoupled), and treat the Higgs masses 
as independent parameters. The three curves vl,v2,v3 
correspond to g?ti induced by the Higgs vertex correc- 
tions with rriA — 100, 200, 300 GeV, respectively, while 
the curve ml corresponds to dxi induced by the A — H 
mixing with m A — 100 GeV (for m A > 150 GeV this 
effect does not impose any considerable constraints). At 
tan j3 > 15 the generated EDM is comparable to the ex- 
perimental limit and at tan/3 ~ 60 it exceeds the ex- 
perimental limit by up to two orders of magnitude. We 
observe that for this choice of the parameters, the sen- 
sitivity to </> M + <f>A t is better than that in Eq.(pl]) by a 
factor of a few. The sensitivity to </> M + ^3 will be as good 
or better if M3 < 60u. We remark that the considered 
parameter space is constrained by the observed B — > S7 
branching ratio. When the effect of the SUSY threshold 
corrections is taken into account, the constraints become 
rather weak 

We further compare the Higgs-mediated EDMs with 
the two-loop effects considered in Ref. [0 . The strongest 
constraint was obtained from the electron or, equiva- 
lently, the thallium atom EDM (the neutron EDM im- 
poses weaker constraints, at least in the naive quark 
model). At low tan /3 the two-loop diagrams are more im- 
portant than the ones considered here, while at tan (3 > 
35 the Higgs induced four- fermion operators provide 
stronger constraints (see Fig. 2 of Ref. B subject to the 
factor of two correction). This statement, of course, de- 
pends on m A and the soft masses. Concerning the n%A 
dependence, the EDM due to the Higgs vertex scales 
as l/m A and the EDM due to the A — H mixing falls 
off faster than l/m\. The Barr-Zee contribution scales 
down roughly linearly with m,A |§ and therefore is more 
significant at heavier ttia- On the other hand, the Higgs- 
mediated EDMs become dominant in the case of heavier 
superpartners. 

In Fig. 4 we present the behavior of different EDM con- 
tributions in this "decoupling" limit, which leads us basi- 
cally to the pattern of SUSY breaking considered before. 
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In order to observe a cross-over from the two-loop contri- 
butions to the Higgs mediation in dxi , we fix tan j3 = 40 
and scale the SUSY mass parameters by the common fac- 
tor X: \A t , b \ = = IX TeV, Arg(A t , 6 M) = tt/2, and 
m sq = 0.6X TeV. 
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FIG. 3. The CP-odd electron-nucleon coupling Cs as 
a function of tan/3. The curves vl,v2,v3 correspond 
to Cs induced by the Higgs vertex corrections with 
m A = 100, 200, 300 GeV. The curve ml corresponds to Cs 
induced by the A — H mixing with rn A = 100 GeV. The SUSY 
parameters are \A t ,b\ = |m| = 1 TeV, Arg(A t ,bl^) = 7r/2, and 
m sq = 0.6 TeV. 
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FIG. 4. Tl EDM contributions as a function of the SUSY 
scale factor X. The curves bl,b2,b3 and vl,v2,v3 correspond 
to the Barr-Zee type and the CP-odd Higgs vertex contribu- 
tions to the EDMs for m A =100,300,500 GeV, respectively. 
The curve ml represents the contribution due to the A — H 
mixing with m A = 100 GeV. The SUSY parameters are 
tan/3 = 40, \A t , b \ = \fj,\ = IX TeV, Arg(A t , b [i) = tt/2, and 
m sq = 0.6V TeV. 

While the Barr-Zee contributions decouple rather 
quickly, the effects we consider stay constant in the limit 
of heavy superpartners. Finally, we note that there are 
the original Barr-Zee contributions to the EDMs with 



the top and bottom quarks in the loop. In SUSY mod- 
els, these are 3-loop effects and do not lead to significant 
constraints (see |18 for the analysis in general 2HDMs). 

To summarize, we have considered the EDMs of heavy 
atoms induced by the Higgs exchange in SUSY models 
with CP violation in the supersymmetric sector. The 
mechanism of Higgs-mediation is insensitive to an overall 
scale of the superpartner masses and grows as tan 3 f3 for 
fixed values of ttia.h- It provides significant and free 
of large hadronic and nuclear uncertainties constraints 
on the SUSY CP-phases for tan f3 > 10 and electroweak 
scale rriA h- 
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